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The complete point symmetry group of the barotropic vorticity equation on the sphere is deter- 
mined. The method we use relies on the invariance of megaideals of the maximal Lie invariance 
algebra of a system of differential equations under automorphisms generated by the associated point 
symmetry group. A convenient set of megaideals is found for the maximal Lie invariance algebra of 
the spherical vorticity equation. We prove that there are only two independent (up to composition 



^vq ■ with continuous point symmetry transformations) discrete symmetries for this equation. 
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1 Introduction 



Lie symmetries of the inviscid barotropic vorticity equation on the rotating sphere (sBVE) were 
computed in [2] and used in [21 [5] in order to derive a point transformation mapping the equation 
in a rotating reference frame to the equation in the rest frame and to construct exact solutions. 
Therein it was also indicated that besides an infinite-dimensional maximal Lie invariance group 
the barotropic vorticity equation on the sphere admits two independent (up to composition with 
each other and with continuous transformations) discrete symmetries, which merely alternate 
signs of two pairs of variables, (the time, the longitude) and (the latitude, the stream function), 
respectively. However, no systematic derivation of discrete symmetries or, more generally, the 
complete point symmetry group of the sBVE was given in the literature up to now. 

With the present paper, we aim to complete the description of point symmetries admitted by 
\ the sBVE by computing the complete point symmetry group Gq of this equation. We simplify 

the computation within the framework of the direct method via combining it with an advanced 
version of the algebraic approach originally proposed in [8j [9], essentially modified in [3] and 
then developed in [HE]. As a result, we prove that in fact the group Gr> is generated by Lie 
£^ ■ symmetry transformations of the sBVE and the above two discrete transformations. 

CN ■ The sBVE is an appropriate equation to demonstrate advantages of the enhanced version of 

the algebraic approach. In principle, the group Gq, might be computed using merely the direct 
method based on the definition of the set of transformations to be found and the prolongation 
of finite transformations to derivatives by the chain rule. This approach is widely applied in the 
literature for finding complete point symmetry groups of single systems of differential equations 
or equivalence groups and equivalence groupoids (i.e., sets of all admissible point transforma- 
tions) of classes of such systems, see, e.g., [TUJ QT1 HI HI LTHl EU [2D] and references therein. 
At the same time, the sBVE is a third-order nonlinear partial differential equation for a single 
scalar function of three independent variables, and one of the independent variables explicitly 
appears in the equation. Hence the mere application of the direct method for the computation 
of the group Gci is too cumbersome. Moreover, the sBVE admits an infinite-dimensional max- 
imal Lie invariance algebra and therefore the version of the algebraic method proposed in [HI [U] 
is not applicable here, because it relies on the computation of automorphism matrices, which 
properly works only in the finite-dimensional case. From the physical point of view, the sBVE 
is of superior importance as it is capable of describing qualitatively the large-scale behavior of 
the flow in the middle of the troposphere. Due to its relevance for the larger atmospheric scales, 
it is especially convenient to consider the vorticity equation in spherical coordinates. 
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Our paper is organized in the following way: In Section [2] we recall the known point sym- 
metries of the sBVE. In Section [3] we describe a method that can be applied to determine the 
complete point symmetry group of a system of differential equations possessing a nontrivial Lie 
invariance algebra. This method involves the notion of megaideals of Lie algebras and properly 
works even for systems whose maximal Lie invariance algebras are infinite dimensional. Sec- 
tion U] is central. After determining a convenient set of megaideals for the maximal Lie invariance 
algebra of the vorticity equation on the sphere, we derive a maximal system of constraints for 
elements of the group Gn, which are related to properties of the adjoint action of Gn on its Lie 
algebra, and then complete the computation of Gn by the direct method using the constraints 
derived within the framework of the algebraic approach. We briefly sum up our results in the 
conclusion. 



2 The model 

Introducing the stream function in the incompressible two-dimensional (i.e. horizontal) Euler 
equations in spherical geometry leads to a single third-order nonlinear partial differential equa- 
tion for the stream function, which is referred to as the barotropic vorticity equation on the 
sphere. It reads 

Ct + ^aCm - VvCa + 2Jty A = 0, C := ^A + ((1 - C 1 ) 

where A and if are the longitude and latitude, respectively, and /i = sin</?, ip is the stream 
function generating an incompressible two-dimensional flow on the sphere, which is related to 
the vorticity £ by means of the Laplacian on the sphere and f2 is the constant angular velocity of 
the rotating sphere. The derived latitudinal variable fi runs from —1 (South Pole) to 1 (North 
Pole). For convenience, we assume the mean radius of the Earth to be scaled to one. 

It was shown in [21 [5] that Eq. ([1]) admits the infinite-dimensional maximal Lie invariance 
algebra, which is denoted as Sq. This algebra is generated by the vector fields 

v = td t - (ip - - ntd x , v = d u z{ g ) = g {t)d^ Ji = d x , 

Ji = /* — 7= » ^A + — ((1 - n^dy. + ndip) , 

„ cos(A + £lt) n sin (A + fit) 9 . 

j 3 = » \ > d x V ' ((i - /i 2 )^ + nd f ) , 

V 1 "A* V 1 - M 

where the parameter- function g traverses the set of smooth functions of t. The structure of the 
algebra Sq is so(3) © (02 C (Z(g))), where the three-dimensional orthogonal algebra so(3) is 
realized by the vector fields J7i, i = 1,2,3, 32 = ^PiV) is a realization of the two-dimensional 
non-Abelian algebra and {Z{g)) is an infinite-dimensional Abelian ideal in Sq. 

An important property of the above family of Lie invariance algebras parameterized by the 
angular velocity fi is that it is not singular with respect to the parameter fi at f2 = 0, i.e. 
it includes the case of the rest reference frame, and it is natural to denote the maximal Lie 
invariance algebra of Eq. ([1]) with fi = by Sq°. Indeed, it was also shown in [21 [5] that Q can 
be set to zero in the above algebra by means of the transformation 

t = t, fl = fx, A = A + £lt, tp = ip — fi/i. (2) 

The same transformation also allows one to set fi to zero in the vorticity equation ([1]) . Note that 
the transformation ([2]) was originally derived in [T3] , where it was used to transform the vorticity 
equation into a reference frame with vanishing angular momentum. The transformation ([2]) can 
also be found by noting that the algebras Sq and are isomorphic and by constructing the 
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mapping relating these two algebras. For our purpose this transformation is especially convenient 
as it leads to a simplified form of both the vorticity equation ([1]) and the maximal Lie invariance 
algebra Sq (i.e. we can work with S^°). Moreover, setting £1 to zero makes clear the physical 
meaning of basis elements of the algebra Sq. If £1 = 0, the vector field D generates simultaneous 
scalings in t and tp and the vector field V is associated with time translations. The elements 
of the form Z{g) are the infinitesimal counterparts of gauging of the stream function up to a 
summand being a smooth function of t. The vector fields Ji, i = 1,2,3, generate rotations 
represented in angular coordinates. 

It was claimed in [21 E] that among continuous symmetries generated by elements from the 
maximal Lie invariance algebra (or, equivalently, <Sg°), there are two discrete symmetries 
admitted by Eq. (fT]) and these are given by the changes of the signs, (t, A, fi, tp) \-t (— t, —A, fi, ip) 
and (t,X, fjb,if}) I—)- (t,X,—fjb,—ip), respectively. While it is straightforward to check by direct 
substitution that these transformations are indeed symmetries of $1]), it is more elaborate to 
derive them directly from the invariance criterion. This was not done in [21 [5]. It is even harder 
to prove that there are no other independent (up to composition with each other and with 
continuous symmetry transformations) discrete symmetries than these two mirror symmetries. 
It is the purpose of this paper to show by determining the complete point symmetry group of 
the sBVE that there are indeed only these two discrete symmetries. 

3 How to find the complete point symmetry group 
via the algebraic method 

It is considerably more difficult to find discrete point symmetries of a system of differential 
equations than its Lie symmetries. The reason for this is that the powerful infinitesimal sym- 
metry criterion is only applicable for transformations depending on continuous group parame- 
ters [Ul El Q3]. This is also the reason why to date no existing computer algebra package, such 
as [3 QjJl [22] can be used for this purpose, because all such packages rely on the integration of 
the infinitesimal determining equations, which by definition only exist for Lie symmetries and 
are linear. Finding discrete point symmetries or the complete point symmetry group of a system 
of differential equations therefore has in fact to be done by hand, using computer programs only 
for related routine calculations in dialog mode. 

Before we start the computation of the complete point symmetry group Go for Eq. (TjQ) with 
fl = 0, let us recall about the general method proposed in [3], which in some sense can be seen 
as a refinement of the technique suggested in [HI [9] by involving the notion of megaideals [T7] . 

Namely, we use the following property: Given a system of differential equations C, for any 
transformation T from the maximal "point symmetry (pseudo) group G of the system C the linear 
mapping T* : — > Q generated by T on the maximal Lie invariance algebra g of the system C 
via push-forwarding of vector fields in the space of system variables is an automorphism of g, 
% € Aut(g), and hence it preserves all megaideals of q. 

The correspondence T — > T* defines a representation of G on g, which is often unfaithful. 
In particular, this is the case if the group G (resp. the algebra g) has a nontrivial center. The 
representation image G* is a subgroup of the automorphism group Aut(g), which may be smaller 
than the entire group Aut(g). The continuous point symmetries of the system C generate, via 
finite compositions, a connected normal subgroup U of G and induce mappings on g which 
can be considered as internal automorphisms of g and which generate a normal subgroup of 
Aut(g). We denote this subgroup by Int(g). Elements of the factor group G/U (more precisely, 
their representatives in G) are interpreted as "discrete symmetries" of the system C which are 
independent up to composing with continuous symmetries of C. As the representations of G 
and U on g via push-forwarding of vector fields are not necessarily faithful, there is no assurance 
on the existence of a bijection between the factor groups G/U and Aut(g)/Int(g). The rigorous 
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consideration gives rise to a number of difficult problems which concern the relation of algebras 
of vector fields and (pseudo) groups of transformations in the infinite-dimensional case and which 
are out of the subject of this paper. At the same time, these problems can be easily solved in 
particular cases, i.e., related to models of fluid dynamics and meteorology. 

If the algebra g is not low dimensional, the computation of Aut(g) itself may be a complicated 
problem. Moreover, the group Aut(g) may be much wider than G*, especially if the algebra g 
is infinite dimensional. If this is the case, in the course of the construction of Aut(g) we will 
spend efforts for finding elements from Aut(g) \ G>, which are in fact needless for determining 
G*. To avoid such needless computations, instead of the condition G> C Aut(g) we can use the 
weaker condition that G*\ C i if t is a megaideal of g. 

In general, a megaideal i of a Lie algebra g is a vector subspace of g that is invariant under 
any mapping from the automorphism group Aut(g) of g [U El; i- e -) = z for any z € i and 
any X G Aut(g). Every megaideal of g is an ideal and a characteristic ideal of g. A set of 
megaideals of g can be computed without knowing of Aut(g). Both the improper subalgebras 
of g (the zero subspace and g itself) are (improper) megaideals of g. If ii and \2 are megaideals 
of g then so are ii -H2, h ni2 and [11,12]) i-e., sums, intersections and Lie products of megaideals 
are again megaideals. If i2 is a megaideal of ii and ii is a megaideal of g then i2 is a megaideal 
of g, i.e., megaideals of megaideals are also megaideals. All elements of the derived, upper and 
lower central series of g, including the center and the derivative of g, as well as the radical and 
nil-radical of g are its megaideals. In order to have a sufficient store of megaideals, we need one 
more way for finding new megaideals from known ones. 

Proposition 1. // io, ii and \2 are megaideals of g then the set s of elements from io whose 
commutators with arbitrary elements from ii belong to 12 is also a megaideal of g. 

Proof. It is obvious that s is a linear subspace of g. Consider an element zq G io such that 
[zo,zi] G i2 for arbitrary z\ G ii. Then for arbitrary T 6 Aut(g) and arbitrary z\ G i\ we have 
[Tzoj^i] = [Xzo, TT _1 zi] = 1[zo, r £~ 1 z\] G i2 as 1~ 1 z\ G ii, and hence [zo,1 _1 .2i] G \2- This 
means that Izo G s, i.e., s is a megaideal of g. □ 

As the megaideals ii and i2 are necessarily usual ideals and hence [io,ii] C t n ii, it in fact 
suffices to consider the case when \2 is contained in io H ii- If io fl ii = {0}, the megaideal s 
coincides with io. A particular case of Proposition [1] with ii = g and i2 = {0} implies that the 
centralizer of every megaideal is a megaideal. 

Once a convenient set of megaideals is found, one should try to obtain maximal restrictions on 
the form of a point symmetry transformation T, which are feasible using the algebraic approach. 
If all derived restrictions on T are taken into account, one has to substitute the restricted form of 
a general point transformation into the initial system of differential equations and proceeds the 
computation of the complete point symmetry group within the framework of the direct method. 

4 Computation of the complete point symmetry group 

We fix the value = 0. It is straightforward to compute the following megaideals of g = S^°: 

g' = (J 1 ,J 2 ,J 3 ,V,Z(g)), g" = (J 1 ,J2,J S ,Z(g)), g'" = 

C g (g') = Z g , = <Z(1)>, Cg (g") = V = (2(9)), C 9 (q">) = (V,V, Z(g)), 

(C s ( 9 '")y = (V,Z(g)), 

where a', Z a and C a (b) denote the derivative and the center of a Lie algebra a and the centralizer 
of a subalgebra b in a, respectively. 

Now we apply Proposition [1] to the case io = k = (C g (g'")) / and vary i2. If i2 = {Z(l)) we 
obtain s = {Z(l), Z(t)) and hence this is a megaideal. We reassign the last 5 as i2 and iterate 
this procedure, which results in the series of megaideals (Z(l),Z(t), . . . ,Z(t n )), n G Nq. 
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Megaideals of which are sums of other megaideals are not essential for the computation of 
the complete point symmetry group Go of the vorticity equation ([1]) for Q = by the algebraic 
method since they give weaker constraints for components of point symmetry transformations 
than their summands. Even if a megaideal is not a sum of other megaideals, its usage may 
result to constraints all of which are consequences of constraints derived in the course of the 
consideration of other megaideals. In order to simplify the computation, we choose a minimal 
set of megaideals which allow us to easily obtain a maximal set of constraints available within 
the algebraic framework. We selected such megaideals from the above list: 

(Z(l)), (Z(l),Z(t)), (V,Z(g)}, (J U J 2 ,J3). (3) 

The general form of a point transformation that can be applied to the vorticity equation on 
the sphere (pQ) for $7 = is 

T: (t,A,/2,$ = (T,A,M,tt), 

where T, A, M and are regarded as functions of t, A, \x and tp, whose joint Jacobian J does 
not vanish. To derive a constrained form of T, we use the selected four megaideals © of S^°. 
For the transformation T to be qualified as a point symmetry of the vorticity equation on the 
sphere, its counterpart 71 push-forwarding vector fields should preserve each of these megaideals. 
Moreover, for any megaideal m of g the mapping induced by T on m is an automorphism 
of m. This property is convenient to use for finite-dimensional megaideals. Thus, the megaideal 
(iTi, J2, 1J3) is isomorphic to the algebra so (3), whose automorphism group is exhausted by 
internal automorphisms and hence isomorphic to the special orthogonal group SO(3). 
As a result, we obtain the conditions 

%Z{1) = T^di + A^d- X + M^<9 A + tf^ = cZ(l), 

%Z{t) = t(T^d~ t + A 4 ,d~ x + + tf^) = d x Z(l) + d Z(l), 

%V = T t 8~ t + A t d x + Mtdf, + = a x V + Z{g), 
3 

T*iJi = ^ bij J j , i = 1 , 2, 3 , 
j'=i 

where g is a smooth function of t which is determined, as the constant parameters c, do, d\, a\, 
0,2, 0,3 and b^, by % and the vector field from the corresponding left-hand side, and (i>ij)ij=i,2,3 
is a special orthogonal matrix. 

We will derive constraints on T* by sequentially equating the coefficients of vector fields in 
the conditions (|4ap ~ (|4d|) and by taking into account the constraints obtained in previous steps. 

Thus, condition (j4a| directly implies that = A^, = = and Vf^ = c. Then the last 
value is nonzero since the Jacobian J does not vanish. The equation ct = d\t + do derived from 
condition (|4bp gives that d\ 7^ and hence the i-component of the transformation T depends 
only on t and the dependence is affine, t = T(t) = cd± 1 t — d^d^ 1 - Condition (j4"cj) is split into the 
equations Tj = a\ (and hence a\ = cd^ 1 7^ 0), Aj = Mt = and = g. Collecting coefficients 
of in condition ()4d|) . we obtain that ^\ = *5>^ = 0. The integration and arrangement of all 
the above equations for the components of T results in the representation 

T = ait + a , A = A(A,/x), M = M(A,/x), = cip + f (t) , 

where 01, ao and c are arbitrary constants with a\c 7^ 0, / is an arbitrary smooth function of t, 
the pair of the smooth functions A and M has nonvanishing Jacobian and additionally satisfies 
equations implied by condition (|4dj) . Up to internal automorphisms of the algebra tS^ which are 
generated by the rotation operators J\, J2 and J3, we can set the matrix (6^)^=1,2,3 to be the 
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(4a) 
(4b) 
(4c) 

(4d) 



unit matrix. Then we obtain the following system of equations with respect to the functions A 
and M: 



M M 
JiA = 1, J2A = = == sin A, J 3 A = = — cos A, (5a) 

J 1 M = 0, J 2 M = y/l- AT 2 cos A, J 3 M = -y/l - M 2 sin A, (5b) 

The equations J ± A = 1 and J X M = imply that A = A + T(/i) and M = M(n). We 
substitute these expressions into the last two equations of (|5b|) and split them with respect to A. 
This gives the conditions y/l - M 2 sin T = and y/l - /i 2 M M = \/i - A/ 2 cosT. As / 0, 
we have that sinT = 0, i.e. T = irk, where k G Z. The same procedure applied to the last two 
equations of (f5"aj) results in the condition 

H (-l) k M 



y/T^J/ 2 y/l ~ M 2 ' 



which is equivalent to M = (— l) fc /i. Then the equation y/l — ^ 2 M^ = y/l — M 2 cos T is iden- 
tically satisfied. 

There are no more constraints which can be derived within the framework of the algebraic 
method. The further consideration is based on the direct calculation of transformed derivatives, 
which is quite easy since the expressions for the transformation components have already been 
specified. Thus, the transformed left-hand side of the vorticity equation ([I]) with £1 = 

Ct + (i'xCn - 4>a(~ x ) = —Ct + {-i) k c 2 {^ - VvCa) 

identically vanishes for each solution of ([1]) if and only if c = (— l) k /a\. This means that up to 
rotations, which are generated by vector fields from { J7i , J2, any transformation from the 
group Go takes the form 

t = ait + a , A = A, u = e/x, i> = —ip + f (t) , 

ai 

where ao and a\ are arbitrary constants with a\ 7^ 0, e = ±1 and / is an arbitrary smooth 
function of t. (We neglect the shift of A by irk as it is a rotation associated with J\ and 
denote (— l) k by e.) A transformation of the above form belongs to the connected component 
of the unity in Go if and only if a\ > and e = 1. Therefore, there are only two discrete 
transformations in Go which are independent up to combinations with each other and with 
continuous transformations. These are, e.g., the transformations with (ai,e) = (—1,1) and 
(ai,e) = (1,-1), where in both the cases we set ao = and / = 0, which merely alternate the 
signs of the variables {t,tp} and {//,?/>}, respectively. 

The transformation just alternating the signs of the variables {A, [i\ is a continuous symmetry 
transformation of the vorticity equation ([1]) for £1 = as in fact it is the rotation by the angle tt 
with respect to the axis corresponding to A = and \i = 0. The above symmetry transforma- 
tions alternating signs of different sets of variables can be combined in order to obtain other 
pairs of simple discrete transformations which are independent of each other up to continuous 
transformations. An example of such a pair is given by the transformations merely alternating 
the signs of the variables {t,\} and {n,ip}, respectively. These transformations coincide with 
those stated in [2J[S]. This completes the description of the complete point symmetry group of 
the barotropic vorticity equation on the sphere with = 0. 

By usage of the transformation (|2|) the above discrete transformations can also be transferred 
to discrete symmetries of the vorticity equation on a constantly rotating sphere. 

Summing up the above consideration, we obtain the following assertion. 
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Theorem 1. The complete point symmetry group of the barotropic vorticity equation on the 
sphere (pQ) is generated by one-parameter groups associated with vector fields from the algebra Sq 
and two discrete transformations, e.g., 

(t, A, /x, tp) i— > (— i, — A, fi, ip) and (i, A, /i, VO i->- (t, \,—fj,,—ijj). 

Corollary 1. The factor group of the complete point symmetry group of the barotropic vorticity 
equation on the sphere (Q]) with respect to its connected component of the unity is isomorphic to 
the group 1*2 x Z2. 

5 Conclusion 

In this paper we verified the claim raised in [2j [5] that the barotropic vorticity equation on 
the sphere possesses two independent (up to composition with each other and with continuous 
symmetry transformations) discrete symmetries. The computation involved two parts, an alge- 
braic step and a step related to the direct method of finding point symmetries. In view of the 
structure of the maximal Lie invariance algebra S^°, we were able to find a sufficiently large 
number of megaideals of 5^° and then selected those of them which were essential for our consid- 
eration, i.e. the megaideals ©. This allowed us to derive important restrictions on the form of 
point symmetry transformations and therefore strongly economized the remaining computations 
which were necessary to be carried out using the direct method. We should in particular stress 
that by taking into account all the constraints that are derivable by the algebraic method, we 
already obtained a strongly restricted form of the admitted point symmetries. Only a single 
constraint, which relates the constants a±, e and c, could not be found from the transformation 
behavior of the megaideals and consequently had to be determined using the direct method. As 
the sBVE is a complicated third-order nonlinear partial differential equation in (1 + 2) variables, 
not deriving the above restricted form would have rendered it quite problematic to compute the 
complete point symmetry group using only the direct method. 

By Proposition [1] we also extended the number of possibilities to determine megaideals of 
Lie algebras. This will be crucial for the computation of the complete point symmetry group 
of other systems of differential equations as the method we proposed in [1] and applied in this 
paper heavily relies on the availability of a large number of megaideals of the associated maximal 
Lie invariance algebras. 

Another novel feature of the present paper is the combining of a simplification of auto- 
morphisms via factoring out internal automorphisms as originally proposed in HH [9] with the 
algebraic technique based on megaideals. This is advantageous for the case under consideration 
as the rotations from SO(3) in angular coordinates have a rather cumbersome representation, 
i.e. already the direct integration of the Lie equations associated with elements of so (3) is a 
nontrivial problem. If the calculation of the complete point symmetry group Go would be done 
without factoring out internal automorphisms, the integration of the Lie equations would be 
implicitly repeated during the computation, which would considerably complicate the calcula- 
tions within the algebraic method. As so (3) is both a direct summand and a megaideal of S^ 3 , 
the extension of any automorphism of so (3) to the complement of so (3) in Sq 2 by identity is an 
automorphism of S^°. Moreover, any such automorphism of is internal as the automorphism 
group of so(3) coincides with the group of internal automorphisms. Hence we can easily factor 
out such automorphisms assuming in the course of the computation that the basis elements Ji, 
J2 and J% are identically transformed. Factoring out other internal automorphisms does not 
essentially simplify the consideration. 

To conclude, it often happens that some discrete symmetries of a system of differential equa- 
tions are known but it is difficult to prove that there are no other discrete symmetries. It will 
therefore be instructive to test the refined algebraic method for the computation of discrete 
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symmetries as presented in this paper with equations which are known to possess nontrivial 
discrete symmetries, such as the potential fast diffusion equation vt = v xx /v x , cf. j!5j . 
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